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ABSTRACT 

We study the prospects for detecting neutrino masses from the galaxy angular power 
spectrum in photometric redshift shells of the Dark Energy Survey (DES) over a vol- 
ume of ~ 20 h~ 3 Gpc 3 , combined with the Cosmic Microwave Background (CMB) 
angular fluctuations expected to be measured from the Planck satellite. We find that 
for a A-CDM concordance model with 7 free parameters in addition to a fiducial neu- 
trino mass of M v = 0.24 eV, we recover from DES&Planck the correct value with 
uncertainty of ±0.12 eV (95 % CL), assuming perfect knowledge of the galaxy bias- 
ing. If the fiducial total mass is close to zero, then the upper limit is 0.11 eV(95 % 
CL). This upper limit from DES&Planck is over 3 times tighter than using Planck 
alone, as DES breaks the parameter degeneracies in a CMB-only analysis. The anal- 
ysis utlilizes spherical harmonics up to 300, averaged in bin of 10 to mimic the DES 
sky coverage. The results are similar if we supplement DES bands (grizY) with the 
VISTA Hemisphere Survey (VHS) near infrared band (JHK). The result is robust to 
uncertainties in non-linear fluctuations and redshift distortions. However, the result is 
sensitive to the assumed galaxy biasing schemes and it requires accurate prior knowl- 
edge of the biasing. To summarize, if the total neutrino mass in nature greater than 
O.leV, we should be able to detect it with DES&Planck, a result with great importance 
to fundamental Physics. 
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1 INTRODUCTION 

Neutrinos are so far the only dark matter candidates that 
we actually know exist. It is now established from so- 
lar, atmospheric, reactor and accelerator neutrino exper- 
iments that neutrinos have non-zero mass, but their ab- 
solute masses are still unknown. Cosmology could provide 
an up per limit on the sum of neutrino masses (for review 
see e.g. lElgarpy fc Lahavl2005l ; lLesgourgues fc Pa stor 2006). 
The growth of Fourier modes with comoving wavenum- 
ber k > k nr will be suppressed because of neutrino free- 
streaming, where 

fc nr = 0.026 ^^y /2 nK 2 h Mpc" 1 , (1) 
E-mail: lahav@star.ucl.ac.uk 



for three equal-mass neutrinos, each with mass m„. 

The current mass upper limit, obtained using Cosmic 
Microwave Background (CMB) WMAP5 data, SN la and 
the BAO from 2dFGRS and SPSS is M„E^m„< 0.61 
eV at 95% CL (|Komatsu et al.ll2009h . The challenge now is 
to bring down reliably the upper limits to the 0.1 eV level or 
even detect the neutrino mass. In this way Cosmology could 
resolve the mass scale of neutrinos. The new generation of 
deep wide surveys can play a key role in setting a tight upper 
limit on the neutrino mass, and possibly detect it if the true 
neutrino mass is sufficiently high. This is due to the order 
of magnitude increase in volume of the new surveys. 

Here we study specifically the ability to set an upper 
limit on the neutrino mass from the galaxy clustering ex- 
pected in the photometric redshift survey Dark Energy Sur- 
vey (DES), combined with Planck CMB measurements. The 
reason a survey like DES would be effective is its large vol- 
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ume, ~ 20 h Gpc , and large number of galaxies, ~ 300 
million. 

Crudely, for a spectroscopic survey, where accurate red- 
shifts are known, the error on the power spectrum scales 
with the survey effective volume V e g as 



2 THE GALAXY SURVEYS 

2.1 The Dark Energy Survey (DES) 



AP(k)/P(k) oc l/\ 



(2) 



On the other hand, the supp ression is proportional in the lin- 

(|Hu et al)ll998l ; iKiakotou et ail 



ear r egime to /, 
2008), where 



/ 



0„ 



93.14ft 2 eV 



(3) 



We expect therefore (when all other cosmological param- 
eters are fixed) that the determination of the upper limit 
on the neutrino mass would be inversely proportional to 
V^eff- From the 2dF Galaxy spectroscopic redshift survey, 
covering a volume of roug hly 0.2 (ft^Gpc) 3 , the upper limit 
on the sum of neut rino mass is about 2 eV at 95 % CL 
l|Elgar0v et al.ll2002f ). Had DES been a spectroscopic survey 
with volume of about 20 (h~ 1 Gpc) 3 , i.e. 100 times larger, 
we would expect an upper limit of 0.2 eV on the sum of 
neutrino mass. Our detailed calculation below yields an up- 
per limit of O.leV for DES&Planck. This is tighter than the 
above back-of-the-envelope calculation, probably as Planck 
priors are incorporated, and the effective volumes above are 
only given crudely. 

However, DES is photometric redshift survey, so the 
radial component of distance to galaxies is significantly de- 
graded, re sulting in a poorer e stimate of the power spec- 
trum (e.g. lBlake fc Bridlell2005l ) . Therefore we prefer in this 
analysis to quantify the galaxy clustering as angular (spher- 
ical harmonic) Ci power spectrum derived in photometric 
redshift shells which are wide enough relative to the pho- 
tometric redshift errors, and to derive the resulting up- 
per limits more carefully and quantitatively. We defer the 
comparison of P(k) and Ce approaches to future studies. 
The utility of photometric redshifts is now well-established, 
with many successful techniques be i ng e mployed (e.g. 
ICollister fc Lahavl 12004 lAbdalla et ail I2009D . The cosmo- 
logical parameter constraints resulting from future photo- 
metric redshift imag ing surveys have been simulated by 
several authors (e.g. | Seo fc Eisensteinl 120031 ; iDolnev et al.l 
12004 I Zhan et "all I2006I ). Applicati on to data such a s the 
SPSS LRG samples were g iven bv lBlake etail [|2007h and 
iPadmanabhan et al.l (|2007 ). These studies mainly empha- 
sized the detection of baryon acoustic oscillations in the 
galaxy clustering pattern. Apart from the specific applica- 
tion to DES, the present paper illustrates more generally the 
determination of neutrino mass from photo-z surveys, to our 
knowledge for the first time. 

This paper is organized as follows. In Sections[5]we sum- 
marize the DES and VHS surveys and Planck, in Section [3] 
we present the photometric redshifts for DES and DES& the 
Vista Hemisphere Survey (VHS) combined filters. Sectionsf?] 
give the formalism for the galaxy angular power spectrum, 
and the associated joint likelihood with Planck. Section [5] 
presents the results for the basic observational and theo- 
retical scenarios, while Section [6] provides extensions of the 
analysis. An overall discussion is given in Section [7] 



The Dark Energy Survey (www.darkenergysurvey.org) is a 
ground-based photometric survey that will image 5000deg 2 
of the South Galactic Cap in the optical griz bands as well 
as the F-band. The survey will be carried out using the 
Blanco 4-m telescope at the Cerro Tololo Inter-American 
Observatory (CTIO) in Chile. The main objectives of the 
survey are to extract information on the nature and den- 
sity of dark energy and dark matter using galaxy clusters, 
galaxy power spectrum measurements, weak lensing studies 
and a supernova survey. This will be achieved by measuring 
redshifts of some 300 million galaxies in the redshift range 
< z < 2, tens of thousands of clusters in the redshift range 
< z < 1.1 and about 2000 Type la supernovae out to red- 
shift z « 1. Observations will be carried out over 525 nights 
spread over five years between 2011 and 2016. The DES vol- 
ume is estimated to be 23.7ft -3 Gpc 3 in the range < z < 2 
assuming a 10ct AB magnitude r < 24 (jBanerii et al.ll2008l ). 

The DES survey area overlaps with that of several other 
important current and future surveys for example the south- 
ern equatorial strip of the Sloan Digital Sky Survey and the 
South Pole Telescope SZE cluster survey. 



2.2 The Vista Hemisphere Survey (VHS) 

The entire DES region will also be imaged in the near infra- 
red bands on two public surveys being conducted on the Vis- 
ible and Infra-Red Survey Telescope for Astronomy (VISTA) 
at ESO's Cerro Paranal Observatory in Chile. In particu- 
lar, the Vista Hemisphere Survey (VHS) (www. vista. ac .uk ) 
will image the entire southern sky (~20000deg 2 ) in the near 
infra-red YJHK S bands when combined with other public 
surveys. About 40% of the total VHS time has been dedi- 
cated to VHS-DES, a 4500deg 2 survey being carried out in 
the DES region over 125 nights in order to complement the 
DES optical data with near infra-red data. The initial pro- 
posal is for the survey to image in the JHK S bands with 
120s exposure times in each band reaching 10cr magnitude 
limits of J = 20.4,// = 20.0 and K s = 19.4. A second pass 
may then be obtained with 240s exposures in each of the 
three NIR filters in order to reach the full-depth required by 
DES. The VHS-DES survey assumes that V-band photom- 
etry will come from the Dark Energy Survey. 



2.3 Planck 

The recently launched ESA's Planck satellite will map the 
CMB with better resolution, sensitivity and frequency (in 
nine bands from one centimeter to one third of a millimeter) 
than previous CMB experiments Q. 

We have obtained forecasts from the Planck satellite 
by u sing the technique described in (|Abdalla fc Rawlingsl 
2007). We have assumed conservative values for the sen- 
sitivity by taking only one usable science channel with 8 
detectors, a noise effective temperature of 120jj,K^/s, an an- 
gular resolution of 10 arcminutes, and 65% sky coverage in 
a one year survey. The input fiducial masses are taken as 



|http : //www . esa. int/ SPECIALS/Planck 
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n b 


h 


C8 


n s 


T 


N u 


0.005 
0.001 


0.255 
0.259 


0.04 
0.04 


0.72 
0.72 


0.9 
0.9 


1 
1 


0.166 
0.166 


2.0 
2.0 



Table 1. The two fiducial models in the Planck simulations. All these 8 parameters vary in the Planck MCMC chain. A flat universe is 
assumed. N v is the number of massive neutrinos. We note the value of the optical depth r assumed here is based on WMAP1, higher 
than the most recent value estimated from WMAP5. 




< \at, m \ >= Ce . 



(4) 



1.0 1.5 2.0 

True Redshift 



Figure 1. Gaussian fits to the true redshift distribution of in each 
of seven photometric redshift slices analyzed in this study. The 
best-fitting Gaussian parameters fj, and a are indicated, where 
p(z) ocexpl-IO-^/Sa 2 ]}. 



Q. v — 0.001 and 0.005, and the other cosmological parame- 
ters considered as listed in Table [T] 



3 PHOTOMETRIC REDSHIFTS 

We determined photo metric redshift estimates as described 
in lBanerii et al. I (|2008l ) using mock sampl es of DES and VHS 
and t he photo-z software package "ANNz" l|Collister fc Lahavl 
120041 ). In brief, artificial neural networks are applied to pa- 
rameterize a non-linear relation between the galaxy redshift 
and the galaxy multi-band photometry. The neural network 
is "trained" using a set of galaxies with known true redshifts 
by minimizing the sum of the squared differences between 
the photometric and true reds hifts. For comp a rison of ANNz 
with 5 other photoz codes see lAbdalla et all (|2009l ). 

Figure [1] shows the distribution of true redshifts in the 
evaluation set for galaxies binned in seven photometric- 
redshift slices, as given in Table [2] We find that these crror 
distributions are well fitted by Gaussian distributions; the 
best-fitting values of the mean and standard deviation for 
the different slices are indicated in the Table. These Gaus- 
sian functions are taken as our model for the redshift distri- 
bution of galaxies in each photo-z slice when analyzing the 
galaxy clustering results. 



4 THE ANGULAR POWER SPECTRUM IN 
PHOTO-Z SHELLS 

4.1 Spherical Harmonic Formalism 

It is common to expand the distribution of galaxies in spher- 
ical harmonic coefficients a^ m , which are then averaged to 
form the angular power spectrum Ce 



(see e.g.|Peebleslll973l ; IScharf et al.lll992l ; IWright et al1ll994l ; 



IWandelt et alj|2001 



The angular power spectrum Ce is a projection of the 
spatial power spectrum of fluctuations at different redshifts 
z, P(k,z), where A; is a co-moving wavenumber. When the 
sky coverage is incomplete the prediction for observed angu- 
lar power spectrum can be estimated through convolution, 
as explained later in the paper. 

We follow here the notation of (|Blake et al.ll2007l '). The 
equation for the projection is: 



C e = — I P (k) ge{k) z dk 

7T 



(5) 



where we assume P(k,z) = Po(k) D(z) 2 , with D(z) the lin- 
ear growth factor at redshift z, and ge(k) contains D(z) as 
defined below. 

We note that this decomposition of P(k, z) is strictly 
only valid in linear theory, and its application at smaller 
scales is an approximation. Moreover, in the presence of mas- 
sive neutrinos P(k, z) cannot be decomposed even in linear 
theory. However, we find that over the redshift z < 2 and 
fc-range of interest P(k z) ps s P (k)D(z) 2 to within 0.5% (cf. 
iLesgourgues fc Pasto j|20o5 iKiakotou et alj|2008h . We have 
also assumed in the above equation a scale-independent and 
epoch-independent bias factor b. Later we allow the biasing 
to vary with redshift. 

Figure [2] shows the expected spherical harmonic Ce (ig- 
noring shot noise) for three values of assumed neutrino mass, 
illustrating that the suppression effect can be measured. Fig- 
ure [3] shows the expected Ct's for three photo-z shells. 

Similarly, C\' 3 is the cross angular power spectrum be- 
tween the two redshift slices i and j , with a kernel ge(k) for 
a redshift slice i: 



Cl' J = 



2 b t b 1 



P {k)g\{k)gl{k)dk 



(6) 



where bi and bj are the linear bias factors for the slices. The 



kernel ge(k) is given by: 
9i{k) = j ji{u) f(u/k)du. 



(7) 



Here, je(x) is the spherical Bessel function and fix) depends 
on the radial distribution of the sources as 



f[z(x)]=p{z)D{z) 



(dx\ _1 
\~dz~) 



(8) 



where x(z) is the co-moving radial co-ordinate at redshift 
z, and p(z) is the redshift probability distribution of the 
sources, normalized such that J p(z)dz — 1. A good ap- 
proximation for equation [5] which is valid for moderately 
large ^>30) is: 
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photo-z 


shell 


galaxy 


DES 


DES&VHS 






fraction 


mean 


(7 


mean 


a 


0.3< z ph 


<0.5 


0.211 


0.405 


0.125 


0.412 


0.145 


0.5< z ph 


<0.7 


0.337 


0.582 


0.125 


0.598 


0.129 


0.7< z ph 


<0.9 


0.215 


0.789 


0.123 


0.805 


0.113 


0.9< z ph 


<1.1 


0.128 


0.975 


0.125 


0.984 


0.116 


l.K z ph 


<1.3 


0.098 


1.203 


0.220 


1.193 


0.142 


1.3< Zph 


<1.5 


0.081 


1.393 


0.260 


1.393 


0.147 


1.5< 


<1.7 


0.027 


1.673 


0.291 


1.593 


0.149 



Table 2. The best-fitting Gaussian parameters /i and a, in true redshift z, where p(z) ex. exp {— [(z — fi) 2 /2a 2 ]} for photo-z shells derived 
from DES alone (5 niters) and DES&VHS (8 filters). 
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Figure 2. The model angular power spectrum for the photomet- 
ric redshift slice z = 0.8 with the exact expression in real space 
(equation [5] with the kernel of equation . The curves are for 
three values of fi„ = 0.000, 0.001 and 0.005, with the other cos- 
mological parameters held fixed at the values given in Table [T] 



Ct 



Po(k=£/x) D 2 (z)x(z) 



~ 2 p{zf 



(dx\~ 
\~dz~) 



dz. (9) 



The above formalism neglects redshift distortions, which we 
examine later. 

The cosmological parameters enter the matter power 
spectrum of fluctuations Po(k), the growth factor D(z) and 
the co-moving distance x(z) in the above equations. For the 
redshift distribution of the sources p(z) we used the Gaus- 
sian functions fitted to the training set data (see Figure [IJ. 
We derived mode l spatial power sp ectra using the "CAMB" 
software package (|Lewis et al.l 2000l). For more technical de- 



tails of these calculations see ijBlake et alj|2007h . 



4.2 Spherical Harmonic Likelihood 

In this work we forecast the measurements, without any 
data, so we use the following relation to express the ratio 
of (log) probabilities between a selected model in the pa- 
rameter space B and the assumed fiducial model A. The 
model describes all the power spectra produced from ran- 
dom Gaussian fields with a certain degree of cross correla- 
tion l|Bucher et al.ll2002l ; lAb~dalla fc Rawlingsll2007n . 



Figure 3. The model angular power spectrum for three photo- 
metric redshift slices centered at z = 0.4,0.8 and 1.6. The curves 
are for f2„ = 0.005, with the other cosmological parameters held 
fixed at the values given in the first row Table [T] The amplitude 
difference between these curves is partially due to the growth of 
structure with redshift and partially due to the variation in the 
photo-z error with redshift. 



log(^ 

PA 



log 



p(aim\B) 
p(ai m \A) 



^^(2Z + 1) (Tr(I- MaMb V )+ l Jm 



log (DetiMAMg 1 ))) 



where / is the identity matrix and the matrices Ma and Mb 
are given by the values of the individual power spectra and 
their cross correlations at a given mode I. As in our case we 
have 7 shells each of them is a 7 x 7 matrix. The diagonal 
terms are 



diag(Mi) = Ci+af 



(11) 



where erf = 



where N / Afl is the average source density 
and /sky = Af2/(47r) denotes the fraction of sky covered by 
the survey. There is no shot-noise in the off-diagonal terms. 

To get insight into this likelihood ratio we note the 
special case of eq. [10] for a shell's auto-correlation (cf. 
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iFisher et al.1l 19941 ): 

\PA ) \P(ai m \A) J 

5 ANALYSIS AND RESULTS 

The joint likelihood is implemented as follows. The Planck 
MCMC chain in 8-dimensional space is centred at one of 
the two fiducial models given in Table [T] The DES likeli- 
hood is evaluated using eq. 1101 with label A correspond- 
ing to the fiducial model, and label B to another point of 
the Planck chain. By evaluating the DES likelihood at the 
Planck chain points we effectively get the joint likelihood of 
DES & Planck. If we are interested in a subset of the pa- 
rameters of interest we can easily marginalize over the other 
parameters by projecting the distribution of the points in 
the resulting chain. 

We first consider a baseline model so we can see "the 
wood for the trees" , where we assume a linear matter power 
spectrum (derived from CAMB), a linear biasing parameter 

6 = 1 for all 7 photo-z slices, no redshift distortion, and we 
ignore the mask. In the discussion below we refer to the 95% 
CL values. 

The first entry in Table is for Planck alone, while the 
second entry is for Planck & DES, both for I — 1,2..., 100. 
We see that the addition of DES improves the upper limit to 
0.11 eV, more than 3 times tighter than Planck alone. The 
corresponding error bar on the high-mass values is 0.13 eV 
from Planck & DES. The likelihood contours in Figures r4l5l6l 
correspond to entry 5 in Table [3] where £ ma x = 300, where 
A£ = 10, which roughly represents the DES incomplete sky, 
with fiducial fi„ = 0.001 and fl„ = 0.005. They show pairs 
of variables (Q u ,as), (fi„,/i), and (f2„,fib), and in Figure 
[7] the probability for fi„, after marginalizing over all other 
seven parameters. Numerical results are given in Table [3] as 
68 % and 95 % upper limits. We note the similarity of results 
of entries 2 and 5 in that Table. We also note that in the 
case of entry 5, the 5-cr upper limit is 0.121 eV, compared 
with the 2-cr value of 0.106 eV. 




2 4 6 8 10 




2 4 6 8 10 12 14 16 



" -3 

v x 10 

Figure 4. The likelihood contours show Q u and <rg from Planck 
alone and Planck & DES , with fiducial f2„ = 0.001 (top) and 
Q u = 0.005 (bottom). The likelihood functions are evaluated for 
(■max = 300, where A£ = 10 (corresponding to entry 5 in Table 
\3\ These plots are for seven cross-correlated DES shells. Contours 
show 68% CL and 95% CL. The sensitivity to the assumed galaxy 
biasing scheme is illustrated with the point O corresponds to b = 
1.00, X corresponds to b = 1.02 constant for all shells, and Y 
corresponds to epoch-dependent b(z) from eg I19l with bo = 1.02. 

coverage by averaging the angular power spectra in multi- 
pole bands of width Ai — 10: 

< cr >= J2w + Wi'/Y^v* + 1} (14) 

£' V 

Entry (4) In Table [3] shows that the upper limit on the neu- 
trino mass is twice as large as for a 'whole sky DES' (entry 2 
in Table [3}. However, if we use harmonics up to t m&x = 300 
(entry 5), then the upper limits are very similar to those in 
entry (2), where ^ ma x = 100. 



6 EXTENSIONS OF THE ANALYSIS 
6.1 Incomplete sky coverage 

When the sky coverage is incomplete the prediction for ob- 
served angular power spectrum can be estimated through 
convolution: 

<C £ obs >=J2 R *,t' c e' (13) 

i 

where the "mixing matrix" Rt^i can be determined from 
the angular power spectrum of the survey window function. 
For a survey of 5000 sq deg or so the mixing matrix smears 
out harmonics ove r A£ m 10 , subject to survey's geometry 
(Blake ct al. 2007|). We illustrate the effect of incomplete sky 



6.2 DES & VHS photoz 

Entries (7) and (8) in Table [3] illustrate the impact of adding 
VHS near infrared photometry, i.e. having photo-z based on 
8 filters. It turns out that the results on neutrino mass look 
very similar to DES alone. This is because the improvement 
of VHS is at high redshift (jBanerii et al.l[2008l ) and the high 
redshift shells have higher shot noise and hence are weighted 
less in the likelihood analysis. 

6.3 Redshift distortion 

These distortions significantly affect the amplitude of the 
projected power spectrum on large scales ^<50, owing to 
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(1) Planck 


1 


100 


1 


M v 


<0.175 


M v < 0.386 


0.034< Mu 


<0.391 


0.004K Mu <0.616 


(2) Planck&DES 


1 


100 


1 


M v 


< 0.067 


M u < 0.112 


0.172< Mu 


<0.304 


0.113< Mu <0.380 


(3) Planck&DES 


20 


120 


1 


M v 


< 0.063 


M v < 0.101 


0.176< Mu 


<0.296 


0.124< Mu <0.362 


(4) Planck&DES 


1 


100 


10 


M v 


< 0.106 


Mu < 0.216 


0.109< Mu 


<0.354 


0.014< Mu <0.476 


(5) Planck&DES 


1 


300 


10 


M v 


< 0.063 


Mu <0.106 


0.169< Mu 


<0.291 


0.113< Mu <0.354 


(6) Planck&DES 


1 


300 


1 


M u 


<0.052 


Mu <0.079 


0.18K Mu 


<0.275 


0.136< Mu <0.322 


(7) Planck&DES&VHS 


1 


100 


1 


M v 


< 0.066 


Mu < 0.110 


0.174< Mu 


<0.304 


0.117< Mu <0.379 


(8) Planck&DES&VHS 


1 


300 


10 


M v 


<0.061 


Mu <0.101 


0.17K Mu 


<0.291 


0.117< Mu <0.353 



Table 3. Table showing derived results neutrino mass Mu for two input fiducial values M v = 0.048 eV and 0.241 eV. The results 
are given for different combination of mock data sets (DES, VHS, Planck) and different range of spherical harmonics l m in,lmax and 
smoothing scale Al. 
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Figure 5. The likelihood contours show Q v and h from Planck 
alone and Planck & DES , with fiducial f2„ = 0.001 (top) and 
fl„ = 0.005 (bottom). The likelihood functions are evaluated 
for tmax = 300 and Al = 10. These plots are for seven cross- 
correlated DES shells. Contours show 68% CL and 95% CL. 




Figure 6. The likelihood contours show Q v and from Planck 
alone and Planck & DES , with fiducial f2„ = 0.001 (top) and 
Qu = 0.005 (bottom). The likelihood functions are evaluated 
for l m ax = 300 and 8£ = 10. These plots are for seven cross- 
correlated DES shells. Contours show 68% CL and 95% CL. 



the relative narrowness of each redshift slice. The amplitude 
of the redshift-space distortions is controlled by a parameter 
(3(z) ~ f2 m (z) ' 6 /b(z) [f], where the quantities on the right- 
hand side of the equation are evaluated at the centre of 
each redshift slice of our analysis. The effect is to introduce 
an additional term to the kernel of equation [7] such that it 
becomes ge(k) + gf{k) where: 



t For a more accurate e xpression in the presence of neutrinos 
see lKiakotou et al.l J2008T) . but this correction is negligible for the 
qualitative argument presented here. 



af(*) = £ J A(u)f'(u/k)du (15) 

jFisher et all Il994l ; iPadmanabhan et all 120071 ; iBlake et all 
l2007f) . 

Figure [8] shows that the redshift distortion is important 
for I < 20. In entry 3 of Table [3] we show the effect of 
analysing 20 < I < 120. Excluding the low £ changes only 
little the results on neutrino mass compared with the case 
where \ow-£ are included (entry 2). 
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Figure 7. The probability for the neutrino mass density Q„, with 
fiducial f2„ = 0.001 and f2„ = 0.005, after marginalizing over the 
other 7 free parameters, for the cases of Planck alone and Planck 
& DES. 
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Figure 8. Real-space (bottom) vs redshift space (top) for the 
nearest (z = 0.4) shell assuming Q v = 0.005 and the other cos- 
mological parameters as in Table [TJ The curves are normalized to 
agree at high I. The assumed biasing is b = 1. 



6.4 Non-linear power spectrum 

Typically, non-linearity in the power sp ectrum appears for 
k > fc max « 0.15 ft Mpc -1 , at z = (e.g. ISmith et al.ll2003l ) 
and higher fc max at higher redshifts. We can estimate the 
equivalent maximum multipole i max of the angular power 
spectrum using the scaling relation fc max = l max /r(z) (see 
equation [9} . For example, for the shells with mean redshift 
z ~ 0.4, 0.8 and 1.6 and assuming a flat universe with f2 m = 
0.25 this corresponds to ^ max ~ 160,290 and 510. Roughly 



Figure 9. linear power spectrum (solid line) vs non-linear 
power spectrum according to Eq. 1171 (dashed line) for the nearest 
(z = 0.4) shell, assuming f2„ = 0.005 and the other cosmological 
parameters as in Table [JJ The curves are normalized to agree at 
I = 1 . The envelope of error bars is derived from Eq. 1181 



speaking it is safe to assume linear theory up to those l? max 
values. 

To check the justification of this simple approximation, 
we considered a non-liner power spectrum in the presence of 
massive neutrinos. The total matter density fluctuation can 
be written as 



5m — fcb&cb + fvSu 



(16) 



where f c b = 1 — fu is the fractional contribution of the 
CDM plus baryon of the present epoch mass density. The 
resulting non-linear power spectrum can be approximated 
i|Saito et al. 1120081 12009| ; llchiki et al.ll2009D as 

P£{k) = fZP?(k)+(l-f v ) 2 P» L (k)+2Ml-f v )P?Uk) .(17) 

where superscripts L and NL denote linear and non-linear 
power spectra respectively, and P^ h „ (fc) is a cross power 
spectrum neutrino and CDM plus baryons. 

This approximation is convenient, as the non- linear 
PS, L (k) can be t aken from fits to N-body simulations 
( ijSaito et al. ll200St ) based on |Smith et al.ll2003l ); "halof it 
= 1" in CAMB). Assuming linearity for the power spec- 
tra which involve neutrino mass is reasonable, as due to 
their free streaming, massive neutrinos remain in the lin- 
ear regime, rather than joining the non-linear evolution of 
CDM plus baryons. We also note that the pre-factor /„ is 
small, which justifies further ignoring non-linear neutrino 
perturbations. For a s lightly different approximation see 
l|Hannestad et alj|2006h . 

For comparison we plot in Figure E] t he errors on the 
CVs (e.g. lDodelso"nll2003l ; iBlake et alj|2007t ): 

a(C e ) = V(2/W(2/ + l)(C t + (18) 

which are far larger than the non-linearity effect for £ < 100. 
As Figure [10] shows the that in this regime the effect of non- 
linearity is less than 5 %. 



6.5 Epoch dependent biasing 

Another source of uncertainty in constrain ing neutrino mass 
from galaxy surveys is galaxy biasing (e.g. lElgar0y fc Lahavl 
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Figure 10. The ratio of C;'s with non-linear to linear power 
spectra according to Eq. 1171 including epoch dependence (solid 
line) for 4 shells, assuming H„ = 0.005 and the other cosmological 
parameters as in Table [T] The curves are normalized to agree at 
2 = 1. We see that up to Z = 100 the effect of non-linear power 
spectrum is less than 5 %. 



2005). The bias systematically increases with redshift for two 
reasons: 

(i) In standard models of the evolution of galaxy clus- 
tering, the bias factor of a class of galaxies increases with 
redshift in opposition to the decreasing linear growth factor, 
in order to reproduce the observed ap proximate constancy 
of the small-scale cluste ring length (e.g. lMagliocchetti et al.l 
l200d : lLahav et alj|2002ft . 

(ii) In a flux limited survey galaxies in more distant 
redshift slices are preferentially more luminous (owing to 
the fixed apparent magnitude thresho ld) and hence more 
strongly clustered (jNorberg et alj|20020 . 

Here we test what happens if the biasing varies with 
epoch . One simple model is the 'galaxy evolving model' (|Frvl 
1996), where the bias evolves as: 

b{z) = 1+ (6o - 1)/D(z) (19) 

where bo is the present galaxy bias for a particular galaxy 
type, and D(z) is the linear theory growth rate. 

The sensitivity to the assumed galaxy biasing scheme is 
illustrated in Figure|4] (bottom). The point O corresponds to 
b = TOO in the fiducial model for all shells, X corresponds to 
b — 1.02, and Y corresponds to epoch-dependent b(z) from 
Eq. [I9]with bo = 1.02. According to that plot, if 6 deviates 
from unity by more than 2% then as changes by 2% and 
Q. v changes by 7%. Therefore biasing is the most sensitive 
quantity in our analysis. In other words, we need to have 
the bias known to 2 per cent accuracy as a larger bias would 
introduce a best fit value outside the one sigma error bar in 
figure 4. 

We see that our results are sensitive to biasing. Fortu- 
nately there are several independent ways of controlling this 
systematic effect : (i) modelling the biasing via halo model 
or semi-analytic simulations, and marginalising over the bi- 
asing parameters; (ii) estimating the neutrino mass from 
galaxy power-spectra derived for different galaxy types and 
checking for consistency; (iii) estimating the biasing empiri- 
cally from weak-lensing map to be produced e.g. from DES 



itself; (iv) estimating biasing from high order statistic such 
as the bi-spectrum. 



7 CONCLUSIONS 

We study the prospects for detecting neutrino masses from 
the galaxy angular power spectrum in photo-z shells in the 
Dark Energy Survey, combined with CMB fluctuations as 
will be measured by Planck. Although the core science case 
for DES is Dark Energy, we see that DES can provide us 
with other important extra science, such as neutrino mass. 
Our main conclusions are: 

• We forecast for DES&Planck a 2-sigma error of total 
neutrino mass AM„ ~ 0.12 eV. If the true neutrino mass is 
very close to zero, then we can obtain an upper limit of 0.11 
eV (95% CL). 

• This upper limit from DES+Planck is over 3 times 
tighter than using Planck alone, as DES breaks the param- 
eter degeneracies in a CMB-only analysis. 

• The results are sensitive to the assumed galaxy biasing, 
and stand if the galaxy bias in known to within 2 per cent. 
This is feasible given other analyses of the galaxy bias su ch 
as the three point correlation function l|Ross et al. 1120071 ). 

• The results are robust to uncertainties in non-linear 
fluctuations and redshift distortion. 

• The results are similar if we supplement DES bands 
(grizY) with the VISTA Hemisphere Survey (VHS) near in- 
frared band (JHK). 

DES can also be used to extract information on neu- 
trino mass via other techniques, e.g. weak gravitational 
lensing, as considered recently for oth er imaging surveys 
pitching et al.ll2008l : llchiki et al.ll2009T l. We note that the 
level of sensitivity for neutrino mass from DES& Planck is of 
much relevance for comparison with the direct measurement 
of the neutrino mass from laboratory experiments. E.g. the 
KATRIN tritium beta decay experiment. Furthermore, the 
DES & Planck measurements can be combined with labo- 
ratory experimen ts to derive more accurate neutrino masses 
l|Host et alj|2007t) . 
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